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ABSTRACT 

We consider evolution of an unstable acoustic mode interacting with an ensemble of stable 
g-modes. We show that the static multimode solution does not exist. We then find the 
condition for the stability of the statistical equilibrium. 

Performing numerical integration of amplitude equations for a simplified system we find 
that the acoustic mode amplitude exhibits a large irregular variability on the timescale given 
by the inverse of the growth rate. 

The g-mode pairs are excited in significantly wider range of detuning parameters than it 
is implied by the parametric instability criterion applied to the average amplitude. However, 
the number of interacting g-mode pairs is reduced because the pairs differing in the detuning 
parameter by less than their damping rates are synchronized and effectively act as a single 
pair. 

We apply the multimode resonant coupling theory to a realistic stellar model. We choose 
a seismic model of the (5-Scuti star XX Pyxidis. Although for some I = 2 modes we find 
amplitudes of the order of a few milimagnitudes, the typical amplitudes of low-degree modes 
are much higher. Taking into account the rotational splitting results in decrease of amplitudes 
by a factor of few which is not enough to obtain consistency with observations. We conclude 
that in this star and likely in all evolved <5-Scuti stars, the resonant mode coupling cannot 
be the dominant amplitude limiting effect. The nonresonant saturation of the driving effect 
must play the role. 

Stars: oscillations - Stars: variables: 5 Scuti 



1 Introduction 

In our understanding of stellar multiperiodic variability, which is the typical 
form of pulsation in main sequence stars and evolved dwarfs, we rely mostly on 
the linear approximation. This suffices if our aim is calculation of oscillation 
mode frequencies. Taking into account nonadiabatic effects allows to determine 
which modes are linearly unstable. This instability usually happens in some 
ranges of spherical harmonic degrees I and eigenfrequencies. There is fairly 
good agreement between theoretical frequency ranges of unstable modes and 
frequencies present in real stars. However, the big unsolved problem is why most 
of the theoretically unstable modes are not seen in observations. This problem is 
related to the more general problem of amplitude limitation mechanism, which 
needs nonlinear approach. 
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The problem with mode ampUtudcs concerns in particular S Scuti stars. De- 
spite the same excitation mechanism as in the giant classical pulsators, namely 
the K mechanism, their pulsational behavior is very different and more diverse. 

First of all, it was noticed already in seventies (Bregcr 1979) that only about 
30% of stars falling into lower instability strip show detectable variability. This 
number is somewhat higher today but still less then one half of the stars in the 
S Scuti domain arc detected as pulsating (Bregcr 2000a). The large number 
of variables with amplitudes close to the detection limit suggest that many 
(perhaps all) of the stars not detected as variable also pulsate but with very 
small amplitudes. 

There is an important subgroup of high-amplitude d Scuti stars (HADS). 
Their pulsational properties resemble those of Cepheid-type variables, i.e., only 
one or two lowest radial modes are present with photometric amplitudes of 
the order of 0.5 mag. These stars are slowly rotating (t;sini<30 km/s) and 
the oscillation amplitudes are constant or almost constant (see, e.g., Rodriguez 
1999). They arc more evolved than the rest of the class. However, the majority 
of S Scuti stars are low-amplitude pulsators with many independent modes most 
of which are nonradial. These typical S Scuti stars can rotate quite fast, with 
vsini up to 250 km/s, but some of them are also slow rotators (Bregcr 2000a). 
The pulsation amplitudes of the dominant modes are below 0.1 mag, and those 
of the remaining modes are usually less than 10 mmag. 

Many low-amplitude S Sciiti stars show long-time amplitude modulations. 
Breger (2000b) analyzed over 30 years of observational data of evolved star 
4CVn. The amplitudes of most modes in this star change significantly in the 
time-scale of ten years or longer. One mode showed rapid decrease of amplitude 
from 15 mmag in 1974 to 4 mmag in 1976 and 1 mmag in 1977 and later started 
growing. There was also a phase jump between 1976 and 1977. The amount of 
amplitude variability in other modes was of the order of 40% during a decade. 

A different example is the star XX Pyx that was observed by Handler et al. 
(2000). The observations covered much shorter time, only 125 nights, but due 
to continuous observation at eight observatories spread aroimd the world and a 
sophisticated method of nonlinear frequency analysis it was possible to detect 
amplitude and phase modulations of many modes with the time-scales from 20 
to 70 days. 

The possible explanation of the apparent difference between pulsational be- 
havior of giant-type pulsators and typical 6 Scuti stars involves the resonant 
mode coupling. The interaction between one higher frequency wave and two 
lower frequency waves is well known in plasma physics (see, e.g., Davidson 1972). 
A detailed analysis of general solutions is given by Wersinger et al. (1980). In 
the case of stellar pulsations it was first studied by Vandakurov (1965, 1979). 
Dziembowski (1982) developed a mathematical formalism suitable to study this 
kind of interaction in S Scuti stars. Application of this formalism to a model 
of ZAMS 6 Scuti star (Dziembowski and Krolikowska 1985) suggested that the 
resonant mode coupling indeed was an effective mechanism of amplitude limita- 
tion. Dziembowski et al. (1988) studied the effects of rotation and showed that 
the amplitudes could be even smaller than in the nonrotating case. However, all 
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these studies concerned a stationary situation. In fact, often happens that the 
stationary solutions are unstable and the dynamical solutions need to be inves- 
tigated. It was done by Moskalik (1985) in the case of three-mode interaction. 
The realistic many-mode case was studied only in the context of tidal-capture 
binaries (Kumar and Goodman, 1996). In that study, the 1 = 2 mode had been 
previously excited due to tidal interaction and the nonlinear couplings caused 
damping of that mode. It is a different situation than in the case of unstable 
acoustic modes in (5-Scuti stars, and we cannot strictly apply that study to our 
case. 



2 Parametric Resonance 

2.1 The Case of a Single G-Mode Pair 
2.1.1 Amplitude Equations 

The amplitude equations describing resonant interaction between an acoustic 
mode with frequency loq and a pair of gravity modes whose frequencies, uji,uj2, 
satisfy the condition |a;o — — -C ^0,1,2 are derived, e.5.,by Dziembowski 
(1982): 

^ = ^,Ao + i^A,A2e--^^-\ (1) 
dr /croio 



^ = 71.2^1,2 + 1 ^ ^ AoAl.e'^'^^, (2) 
dr 2cri^2-'i,2 



where A are dimensionless complex amplitudes of the interacting modes at the 
surface, 7 denote growth rates, / are moments of inertia, H is the coupling 
coefficient, and the dimensionless time, t, and frequencies, a, are related to the 
real ones, t,u! by the relation T/t=u)/a=y/4TTG<~p>. The frequency mismatch is 
given by Aa = cto — ai — (T2 , and e distinguishes the case of two- mode resonance, 
where the subscripts 1,2 denote the same mode, and the case of three- mode 
resonance, where they denote different modes. In the former case e=l/2 and in 
the latter case e= 1. The acoustic mode growth rate is positive and the g-mode 
growth rates are negative. 

Dziembowski (1982) made an explicit use of adiabatic assumption with ar- 
tificially added growth rates. This yields real coupling coefficient. However, 
following a derivation procedure of Van Hoolst (1994b), wc find the nonadia- 
batic amplitude equations in very similar form with the complex H 

= 7o^o + i-^^l^2e-'^'^^ 
dr 2aoIo 

^ = 7i,2Ar,2 + i^^Ao^^,re'^-. 
dr ^cri^2-'i,2 

If wc put 7i = i/e'" and shift complex phases of g-mode amplitudes by —ia/2 
the above amplitude equations are reduced to Eqs. (1),(2). Therefore, without 
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loss of generality, wc will use them with real and positive H. Furthermore, we 
rescale the coupling coefficient 



H 



and the g-mode amplitudes 



^'1,2 = ^1,2 



in order to obtain simpler form of amplitude equations 

= 70^0 + i-;7^i^2e , 



dr 
d^;,2 
dr 



C 



7i,2^'i,2 + i^^o^2,i*e- 



iAo 



(3) 

(4) 

(5) 
(6) 



There are various long timescales which play roles in the evolution of our 
system. These are the detuning scale, |AfT|~^, the nonadiabatic scales, |7|~^ 
and the coupling scale |CAo|~^. 

If the amplitudes are sufficiently small then only linear terms play role and 
the amplitude of the acoustic mode grows exponentially. Eventually it reaches 
the point where the nonlinear terms in Eqs. (6) become important. Assuming 
sufficiently slow growth of the acoustic mode and using the procedure of Dziem- 
bowski (1982) wc obtain a criterion for the parametric excitation of the pair of 
linearly damped g-modes 



l^ol > -Bcr = 



47172 
(02 



Ag 
71+72 



(7) 



When the acoustic mode amplitude exceeds the critical value the g-mode am- 
plitudes start to grow and finally the nonlinear term in Eq. (5) can stop the 
exponential growth of the acoustic mode 

In further applications it will be more suitable to use real form of amplitude 
equations. By substituting 



(8) 



where Bj are the real and positive amplitudes, and c^j are the real phases of the 
considered modes, we obtain 

dSo 
dr 
dBi,2 



C 

= 70^0 + — BiB2sin$, 
C . 



dr 
d$ 

d7 



= 7i,2-Bi,2- — -BoB2,isin$, 



Aa+- 



B1B2 



Bo 



'B2 Bi 
""'^Bl^B-, 



cos$. 



(9) 
(10) 

(11) 



where $ = Acrr + ^0 — (/>! — (j)2- 
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2.1.2 Equilibrium Solutions 

Here, we repeat the analysis of Dziembowski (1982), who derived the formulae 
for the equilibrium amplitudes and the criteria for the stability of the equilibrium 
solution. 

If we set l.h.s. of Eqs. (9)-(ll) equal to zero we find equilibrium solutions 

B'o = + (12) 

Bl,2 = ^J^^^{l + q^), (13) 
cot^^ = q, (14) 

where q = Acr/7, and 7 = (70 + 71+72)- 

The stationary solutions describe possible physical state only if they are 
stable to small perturbations. Linearizing Eqs. (9)-(ll) around equilibrium 
solution (Eqs. 12 14) we have 





" SBo/Bo ' 




To 


-70 


-70 






" 5Bo/Bo ' 


d 


SB^/Bi 




-71 


71 


-71 


-qiv 




SBi/Bi 


d7 


SB2/B2 




-72 


-72 


72 


-qi2 




6B2/B2 




(5$ 




_ 9(270-7) 


9(271-7) 


g(272 - 7) 






(5$ 



(15) 

The solution of such a linear differential equation is the exponential time de- 
pendence exp(AT) and the value of A is given by the characteristic equation 

= A^ + aiA^ + O2A^ + a3A + a4, (16) 

where 

ai = -27, 02 = 7^(1 + g^)-4(7^(7o7i+ 7072+71 72), 

0-3 = 4707172(1 + 3(7^), 04 = -4(1 + (7^)7071 727- 

The equilibrium solution is stable if the real parts of all eigenvalues A are neg- 
ative. This is given by the Hurwitz criteria 

Wi=ai>0, W2=aia2-a3>0, 

W3=a3W2-ala4>0, W4 = a4W3>0. 

It can be easily shown that only first and third criteria are independent in our 
case. They are 

7 = 7o + 7i+72<0, (17) 

473(l+g2)_(i+3g2)j4^^^^^2(l+3g2)+273(l+g2)_8g2^(7o7i+7o72+7i72)]>0. 

(18) 

The first criterion says that driving of the pnmode cannot be higher than damp- 
ing of the g-modes. The second criterion is a quadratic inequality for q"^. It 
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is fulfilled if , and equivalently |Acr|, is confined to a certain range or if it 
exceeds certain value. 

If 7o ^ l7i,2|) the equilibrium acoustic mode amplitude (Eq. 12) is nearly 
equal the critical amplitude (Eq. 7). Moreover, according to Wersinger et al. 
(1980), the terminal state is then a stable equilibrium or a periodic limit cycle 
with average p-mode amplitude close to the equilibrium value. This means that 
the critical amplitude, given by Eq. (7), may always be treated as an estimate 
of the acoustic mode amplitude. 



2.1.3 Energy Equations 

In dynamical considerations it is often useful to use mode energies instead of 
the amplitudes. The general formula for the energy of an oscillating mode is 



The energies in l/2/oc7o units in our case are 

£0 = UoP: 



-fc'1,2 - X 2 l^l.2| ' 
-'0^0 



1 

2"e 



Tr-^1,2' 



(19) 
(20) 



where we used Eqs. (4), (8) and we set 

0'1,2 ~ o'o/2, 

because only the modes satisiying this condition play a role. 
Now we apply Eqs. (9)-(ll) to Eqs. (19), (20) and we obtain 



— = 270^0 + leC^EQEyE-i sin$, 



dr 
dEi,2 
dr 
d$ 
d^ 



— = A(7 + C 



27i,2£^i,2 - C/^J^^sin$, 
Eq 




(21) 
(22) 

cos$. (23) 



From these equations we see that the time derivative of the total energy in 
oscillations is 



^Etotal 

dr 



27oi;o + 27i£;i 



for the two-mode resonance 



270-Eo + 271 £^1 + 272S2 for the three-mode resonance 



(24) 



The coupling terms are only responsible for the energy transfer between the 
acoustic mode and the gravity modes. Since the growth rates have different 
signs it is possible to have the solution with constant total energy, i.e., the 
equilibrium state, or with the total energy varying around some average value, 
i.e., the limit cycle. 
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2.1.4 Equal Damping Rates 

We noted earlier that the modes in aU important g-modc pairs have approxi- 
mately equal frequencies, cti «C72. It turns out that their damping rates have 
the same property, 7if«72=7g. If we assume exact equality we see from Eq. (22) 
that 

d(i;i-£;,) 

= 273(^1-^2)- 

Since 7g<0, the difference between the two energies decays exponentially. Thus, 
we further assume that the energies, as well as the real amplitudes of the g-modes 
in each pair are equal. The energy of the pair is 

Eg = El + E2 = 2Ei 

and then we simplify Eqs. (21)-(23) 

2^oEo + C^EgSm^, (25) 



dEp 

dr 

d^ 

dr 



2jgEg^C^EgSm'S>, (26) 



Identical equations are obtained in the two- mode resonance case, where e= 1/2 
and Eg = Ei. 



2.2 The Case of Many G-Mode Pairs 

The approach adopted by Dziembowski and Krolikowska (1985) was to con- 
sider an ensemble of g-mode pairs that may be excited through the parametric 
resonance instability and to determine probability distribution for the p-mode 
amplitude at onset of the instability. They assumed that the nonlinear interac- 
tion with first excited pair will halt the growth of the p-mode amplitude. This 
could be partially justified when the constant amplitude solution was stable. In 
the case of the ZAMS star model they considered this was often the case. How- 
ever, in her unpublished PhD thesis Krolikowska showed later that in somewhat 
evolved objects the interaction with most likely excited g-modes cannot halt 
the amplitude growth. Therefore, we have to consider the case of many g-mode 
pairs interacting with the p-mode. 
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2.2.1 Amplitude Equations 

The amplitude equations describing the case of many g-mode pairs are obtained 
as a straightforward generalization of Eqs. (1),(2) 

C\A^ rr 

dr 2ct^_2/i,2 

where N is the total number of pairs taken into account. Using the same 
procedure as in the case of a single pair we obtain real amplitude equations 

^ = 7oBo + E^i3iB^sincl>„ (30) 

dB-' C 

= li2Bi2-^BoBl,sm^,, (31) 

Now, let us try to find multimode stationary solution of the system. Putting 
left-hand sides of Eqs. (30)-(32) equal to zero we obtain 



70^0 = -T^^BiBisin^j, (33) 
liA = Y^iiBosm^j, (34) 

AT 



^Bol^ + ^]cos<^j = Aa. + ^y^BlB^cos^k- (35) 
2 °\Bi BiJ ' ^ 2 1 2 V ; 

Prom Eq. (34) we get 



and 



From Eqs. (35)-(37) we obtain 



1 ^ 

(7^ +^') cot$,- = Aaj + -^Y^^HB^f cot$fc. (38) 



Bi — 

fe=l 
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If we multiply Eq. (34) for the mode 1 by the same equation for the mode 2 we 
get 

Eqs. (38), (39) form a set of 2N equations for N + 2 unknowns, namely Bq, 



E 



7f(B*)^cot$fe 



and N values cot<l?j-. This set has in general no solutions if the number of 
equations is larger than the number of unknowns. This happens for N>2. We 
see that the stationary solution exists only if the number of excited pairs is 
not larger than 2. In further sections we will study time evolution of the mode 
amplitudes B and phases $ by means of numerical integration of Eqs. (30)-(32). 



2.2.2 Statistical Equilibrium 

In the absence of stationary multimode solutions we may still expect that the 
system will reach some kind of statistical equilibrium. In this section we find 
a necessary condition for the existence of such a state. We restrict our consid- 
erations to the case of equal damping rates in each pair, 7^ =73 =7j. As we 
mentioned earlier, this is not far from the relevant realistic situation. Then, sim- 
ilarly as in the case of one g-mode pair, we have Bf =B2=Bj and Eqs. (30)-(32) 
become 

N ^ 

loBo + Y^^B^sin^j, (40) 

IjBj-^BoBjsm^j, (41) 

^ C B^ 

+ y cos $fe - CjBo cos . (42) 

Now, we introduce new variables 



dBo 
dr 

dBj ^ 
dr 
d$j 
dr 



X = Bq, Yj = Bj sin , Zj = Bj cos . 



(43) 
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Then, from Eqs. (40)-(42) we obtain 



N 



X = 2joX + Vx^CkYk^Y^+Zl (44) 



Yj ~ ^3 Yj 



(Aa,+^pCkZk^Y,^ + Zl^ , (45) 



N 



The divergence of the flow tells us whether the volume in the phase space 
{X,Yj,Zj) evolved by following the orbits obtained from Eqs. (44) (46) expands 
or shrinks exponentially. If it expands, the system cannot reach the statistical 
equilibrium. After some algebra we obtain the divergence of the flow 

From this we obtain the necessary condition for the stability of the statistical 
equilibrium, which is 

JV 

i=i 

i.e., the summed damping rate has to be larger than the acoustic mode driving 
rate. This is the generalization of the criterion obtained by Wersinger et al. 
(1980) in the case of a single g-mode pair. 

In the same way as in the case of a single g-mode pair we derive equations 
for the time evolution of energies of the interacting modes 



dEo 
dr 



270^0 + /^$^C'j^j sin (49) 



dEj 
dr 



: 2'yjEj-^/EoCjEj&m^j, (50) 

d$- 1 ^ n- 

dr 2y/Eo 

where Ej is the energy of the j-th pair. From Eqs. (49), (50) we see that the 
time derivative of the total energy in the system is given by 

^-f^+i:^-2(>E„+x:7,i^,). (52) 
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If the damping rates are small compared to the growth rate, the number N 
of actively interacting g-modes necessary to halt the acoustic mode growth is 
large, according to Eq. (48). In such a case the system reaches a state with 
time-dependent and limited amplitudes. The average time derivative of the 
total energy (Eq. 52) is then equal to zero. 

If we assume marginal stability, i.e., the sum of damping rates of all active 
g-mode pairs is just equal to the acoustic mode growth rate (see. Eq 48), then 
the total energy derivative (Eq. 52) may be written as 



d-E^total 

dr 



N 



70^0 + J2^j\= 270 (£^0 - E,) , (53) 



where Ej is the mean energy of active g-mode pairs, weighted with the damping 
rates. In the statistic;al ciquilibrium the time-averaged total energy derivative is 
equal zero, which implies 

{Eo) = {E,}. (54) 

This equation looks like an energy equipartition. In fact, the summed damping 
is higher than the driving (Eq. 48), so the mean g-mode pair energy is lower 
than the acoustic mode energy. 

One could think that Eq. (54) could be used to estimate the acoustic mode 
amplitude in the statistical equilibrium. It would be possible if we assumed that 
active were only pairs which are unstable according to Eq. (7) with the acoustic 
mode amplitude averaged over the long time. However, it has to be checked 
numerically whether such an assumption is justified. 



3 Time-Dependent Solutions 

In this section we analyze time evolution of simple systems. This will help us to 
understand behavior of modes in realistic stellar model. We will see that there 
is a qualitative similarity between simple cases and more realistic ones. 



3.1 Single G-Mode Pair 

The case of a single pair of stable modes was studied in details by Wersinger et 
al. (1980). Here we present two examples of such system. 

If the conditions for the stability of the equilibrium solution are satisfied 
the system reaches this equilibrium. More interesting is the case of unstable 
equilibrium and stable periodic limit cycle. An example is presented in Fig. 1. 

Initially, the amplitudes are low and the acoustic mode amplitude grows 
according to the linear driving rate. When the amplitude becomes sufficiently 
high there is significant transfer of energy to the g-mode pair and ultimately 
the acoustic mode amplitude growth is halted. The evolution proceeds rapidly 
on the coupling timescalc. First, nearly all energy of the p-modc is transfered 
to the g-mode pair. Then the energy is returned back to the p-mode, but not 
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-2000 



5 10 15 20 25 30 

time 

Fig. 1. Time evolution of the system of an unstable (acoustic) mode and a pair of stable 
(gravity) modes. The parameters are 70 = 1, 71,2 = —3, Act = 2, C = 1. Upper panel shows 
the amplitude of the acoustic mode. Middle panel shows energies of the acoustic mode (solid 
line) and the g-mode pair (dashed line). Lower panel shows the time derivative of the total 
energy (Eq. 24). The value of 70 sets the unit of time and frequency. The units of amplitudes 
and energies are set by the value of coupling coefficient. 
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completely. Part of the energy is dissipated. The state of constant amplitude 
modulation is reached only after the excess of energy gained in the initial phase 
is dissipated, i.e., after the time about 2O/70. This excess is an artifact of the 
initial conditions. The mean energy balance (Eq. 52) implies 

(55) 

The period of the limit cycle is about 1/70- It is so, because the system 
spends most of the time in the phase of exponential growth of the acoustic 
mode amplitude. Finally, from the bottom panel we see that the averaged time 
derivative of the total energy is indeed equal zero. 

Fig. 2 presents the case when neither a stable equilibrium nor a stable limit 
cycle exist. Such a situation takes place, in particular, when condition given by 
Eq. (48) with A'' = 1 is not satisfied. The amplitude growth is unlimited. This 
growth is modulated in a similar way as in the periodic limit cycle solution, but 
amplitudes in consecutive cycles arc higher and higher. The reason why the 
situation is unstable is that the g-mode damping is low and the pair does not 
manage to loose enough energy to balance the unstable mode driving. In lower 
panel of Fig. 2 it is clearly seen that the average time derivative of the total 
energy is positive. 

We have to add here, that Eq. (48) is not a sufficient condition for the 
stability. Wersinger et al. (1980) showed examples of solutions with unlimited 
growth despite satisfying this condition for N =1 (a single g-mode pair). We 
found similar cases for N>1. What is needed to stop the amplitude growth is 
some excess of the damping. 



3.2 Many G-Mode Pairs 

In this section we study the time evolution of a simplified system of interacting 
modes in the case when several pairs are necessary to halt the acoustic mode 
growth. We set, as previously, 70 = 1 and Cj = 1 for all A'' = 19 pairs. The 
damping rates are a few times smaller than the growth rate and for all g-modes 
we set 7j = — 0.3. The detuning parameters, Aaj, are equally spaced between 
— 13.824 and 13.824 with the separation 1.536. Energies of all interacting modes 
are shown in Fig. 3 and the time derivative of the total energy is shown in Fig. 4. 



At the beginning the acoustic mode amplitude grows exponentially. When 
the amplitude reaches approximately the critical value given in Eq. (7) the 
energy transfer begins. The lowest critical amplitude among the pairs in our 
system is for the one with Act = and this pair is excited first. However, 
this pair alone cannot halt the p-mode growth since the damping rate is too 
low. Gradually, other pairs are excited, but it takes some time before their 
amplitudes are large enough to affect the p-mode growth. The initial excess 
energy is dissipated by modes in a rather wide range of detuning parameters. 

Ultimately, after the time of about I5/70 the system achieves the statistical 
equilibrium in which only 7 closest-to-resonance pairs survive. Let us note that 
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Fig. 2. The same as in Fig. 1 but with parameters that exclude stable solutions, either 
equilibrium, or limit cycle. The parameters are 70 = 1, 71,2 = —0.4, Act = 2, C = l. 
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Fig. 3. The time evolution of tlic multimode system. The energy of the acoustic mode is 
plotted in the back wall of the 3-D plot. The remaining lines denote energies of g-mode 
pairs. The detuning parameters of the pairs are equally spaced and the difference between 
consecutive Aa's is 1.536. Other parameters are the same for all pairs, i.e.,fj = —0.3, Cj = 1. 
The growth rate of the acoustic mode is 70 = 1 ■ This sets the units of time and frequency. 
The unit of energy is determined by the value of coupling coefficients. Initial energies of the 
g-mode pairs are set randomly around the value of 0.01. Initial acoustic mode energy is 0.3. 

the number implied by the necessary stability condition (Eq. 48) is 4. The rapid 
oscillations seen in this phase occur on the timescale determined by the largest 
detuning parameter of the active pairs and this timescale turns out to be crudely 
equal to the nonadiabatic timescale 7(7^. Surprisingly, the result is similar as in 
the case of a single g-mode pair. This result, which is essential for understanding 
of amplitude limitation mechanism will be explained in Section 4. 

The overall behavior of the acoustic mode is characterized by an irregular 
variability with a large amplitude modulations. In Fig. 4 we see that the total 
energy of oscillations is conserved only on average, which means that there is 
energy exchange with the background. 

Fig. 5 presents the energies of active pairs (upper panel) and the energy 
budget for the acoustic mode (lower panel) in a selected time-interval in the 
statistical equilibrium phase. What is important, the local energy maxima of 
various g-mode pairs occur at various times. This lack of synchronization is 
essential. Synchronized pairs act just as a single pair, as we will show in the 
next subsection, and they are not efficient energy receivers. 

In the lower panel we see that the acoustic mode energy balance is zero 
on average. The coupling causes an appreciable energy exchange between the 
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Fig. 4. The time derivative of the total energy in the system shown in Fig. 3. 



acoustic and the gravity modes. The net energy transfer necessary to keep the 
g-modes ahve is relatively small. 

The rapid and nearly complete energy exchange between the acoustic and 
the gravity modes, manifested by very large range of amplitude variations, is 
similar to that in the case of a single g-mode pair. This is so because each 
individual energy exchange act is typically controlled by a single pair. The 
irregularity in the present case is caused by involvement of other pairs. In 
Fig. 5 we see examples when a new pair takes over the control. This happens 
when the previous pair have low amplitude. Which of the remaining pairs takes 
the control depends on the relative phase. It must be favorable, i.e., |Ao|sin$j 
must be large negative. 

From energy balance (Eq. 52) we obtain 

where the summation is done over the active pairs. This is essentially the same 
result as in the case of a single g-mode pair (see Eq. 55). In our example, the 
sum of the damping rates of the active pairs is about 2 times higher than the 
growth rate, 70. Thus, the average g-mode pair energy is about 2 times lower 
than the average energy of the p-mode. 

It is of interest that most of the active g-mode pairs should not be excited 
according to the parametric instability criterion (Eq. 7) applied to the time- 
averaged p-mode amplitude. However, the formula for the critical amplitude is 
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Fig. 5. Upper panel; energies of active pairs in the statistical equilibrium; the zero-levels 
of the energies are shifted proportionally to Acr. Lower panel; driving of the acoustic mode 
(solid line) given by the first term in Eq. (49) and energy transfer to the g-mode pairs due to 
coupling (dashed line), given by the sum in Eq. (49). 
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valid only when the p-mode complex amplitude, Aq, is assumed constant or to 
vary sufficiently slowly. Here this assumption is not satisfied. One can easily 
sec how the varying Aq modifies the parametric excitation criterion. 

Let us assume ^o^^o'^"'^' ^^'^ is constant. Then, Eq. (2) becomes 



dA^ 
dr 



= 71,2^1,2+1 



H 



20-10-^1,2 



^0^2,1' 



i(ACT-|-77)r 



(57) 



i.e., only detuning parameter is modified, 
becomes 



This means the criterion in Eq. (7) 



uil = u„|>B; = 



47172 
C2 



71 + 72 



(58) 



The effective detuning parameter is now Act + 77. In particular, if 77= — Act the 
criterion becomes much less restrictive for the pairs with detuning parameter 
Act than in the case of ^0 = const. 

In the case presented in Figs. (3) (5) the amplitude of the acoustic mode 
varies in a very complicated way. We may treat these variations as consisting 
of many components oscillating with various frequencies rj. A g-mode pair with 
the detuning parameter Act^ is affected mainly by components with ryw — Act^. 
This implies that pairs with relatively large detuning parameters may be excited 
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even when the average acoustic mode amplitude is below the critical value given 
by Eq. (7). Therefore, we cannot estimate average acoustic mode amplitude in 
a simple way described at the end of Section 2. 



3.3 Pairs with Close Detuning Parameters 

When we conducted similar numerical experiments with two or more g-mode 
pairs having very close detuning parameters, AcTj, we observed a synchroniza- 
tion of phases of these pairs. This effect may be demonstrated analytically with 
the use of Eqs. (40) (42). 

Let us consider two pairs, ji, and j2 with Acr^j =Aaj^ , as well as =7j2 =7 
and Cj^ =Cj2 =C Then, from Eq. (42), we have 

^ *--*-.=CBo^^^^^^^^^=5sin^^i^, (59) 



dr 2 "2 



where 



(J) -I- d) 

5 = Cgosin . (60) 



Upon integration we get 



tan ^^1^ = tan ^lilZplk exp ( / 5dr' ) . (61) 



The integrand is a rapidly varying function, but it has a nonzero mean value. 
This value must be less than zero. To see this let us write Eq. (41) for the two 
pairs in the form 

dlnS,-, , C 

=7-^gosin^,-,,,. (62) 

In the statistical equilibrium we have 

C(Bosin$j-,2)=27<0. (63) 

This holds if both phases ^ are in the (— 7r,0) range when Bo is large. Then, 
the average of the two phses, {^j^+'^j^)/2, is also in this range during large Bq 
phase, implying 

(5)<0. (64) 

In fact, the values of the sine functions of both phases as well as of their average, 
are of the same order, except short periods of time when the phases are close to 
0,±7r. This implies the left-hand side of Eq. (64) to be of the order of 27, like 
in Eq. (63), and the synchronization occurs on the timescale |7|~"'^. 

The phase synchronization has an important consequence for the multimode 
coupling. We may replace in Eqs. (49)-(51) the energies of the pairs with the 
same $, 7, C, and Act by the sum of their energies. Then, the equations for 
the other pairs remain unchanged. This means that the synchronized pairs act 
effectively as a single pair and in this way we may reduce the number of pairs 
and the number of equations to solve. 
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In a more realistic situation, the two pairs do not have exactly equal detuning 
parameters. Then, instead of Eq. (59) wc have 

dr 2 2 2 ^ ' 

When we use Eq. (64), then Eq. (65) becomes 

d Act,, -Act,, „ . 

d7 2 ^ + 27sin^l^-^. (66) 

It's solution is 

tan . 



: 1 _ 1^1 _ 22tanh ( - -22|7|r + t) , (67) 

where 2;= (Actji — Acrj2)/(4|7|), and T is the integration constant whose value 
is determined by the initial conditions. If |2:| -C 1 then for 00 we get 

2. = ^^^^^. (68) 

We see now, that the phases are synchronized if the difference of detuning 
parameters is much less than |7|. In fact, we found in numerical experiments 
that the pairs may be treated as synchronized if the Act difference is just smaller 
than 7. Only for larger difference the pairs are really independent. 

The synchronization has important consequence for the stability condition 
(Eq. 48) in which we should count only independent pairs. 



4 General Properties of Multimode Solutions 

In the previous section we studied systems where all pairs had the same damping 
rates and coupling coefficients. It turns out, however, that the main qualitative 
properties of the multimode solutions are valid also in more realistic systems, 
where the pairs have various damping rates and coupling coefficients. 

Firstly, the pairs are excited in significantly wider range of Act than implied 
by Eq. (7) with the acoustic mode amplitude replaced by it's average value. Sec- 
ondly, the pairs with very close detuning parameters are synchronized. Thirdly, 
the p-mode amplitude is very strongly modulated due to nearly complete en- 
ergy transfer between the interacting modes. Finally, the modulation timescale 
is given by the inverse of the growth rate 70. 

The wide range of detuning parameters of active pairs is the result of rapidly 
varying amplitude of the acoustic mode (see Subsection 3.2.2). 

The pairs with close detuning parameters and equal damping rates and cou- 
pling coefficients are synchronized, as was shown in Subsection 3.3. It turns 
out that the pairs with different damping rates and/or coupling coefficients also 
tend to synchronization if their detuning parameters differ less than the smaller 
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of their damping rates. In such a case, in a statistical equihbrium, we have 

(Bosin^j^) 

(see Eq. 63). The synchronization causes both left-hand sides of above equations 
to be roughly equal. If the right-hand sides differ significantly, it is obvious that 
these equations cannot be satisfied simultaneously. Only the one with lower 
value of \-f/C\ can be satisfied, because the sine functions have limited values. 
The energy transfer to the g-mode pair with higher \j/C\ is inefficient and this 
pair gets damped. This means that we may treat such pairs as one efi^ective 
pair, like in the case of equal damping rates and coupling coefficients. Here, 
this effective pair is just the one with the smallest \'y/C\ ratio. 

The large amplitude modulation is the result of nearly complete energy ex- 
change between the p-mode and the g-mode pair that currently controls the 
interac;tion. Tlic; mode energies averaged over all active pairs and over a long 
time satisfy Eq. (56). However, if we have groups of synchronized pairs, each 
group has to be counted as a single pair. We expect, just as in the simplified 
case considered in the previous section, that the average energy of the g-modes 
is a few times lower than the average p-mode energy. 

The fact that the amplitude modulation timescale is of the order of is 
caused by the synchronization of the g-mode pairs with close dc;tuning param- 
eters and the stability criterion given by Eq. (48). In a realistic situation the 
damping rates of the active g-mode pairs are much smaller than the driving 
rate of the acoustic mode. This means that the number of active pairs is high, 
namely it has to be higher than Nmin=lo/\'yj\ (see Eq. 48), where 'fj is the mean 
value of the damping rate of the active pairs. Moreover, the independent active 
pairs have to be separated in the Act space by more than |7j|. This implies the 
range of detuning parameters occupied by active pairs have to be larger than 
Nmin\jj \ =7o- The maximum Act of the active pairs, which gives the shortest 
modulation timescale, is thus higher than 7o/2. In fact, it should be typically 
equal to a few times 70. 

5 Application to a Realistic Stellar Model 
5.1 The Model 

Numerical study of oscillation in a realistic situation must begin with model 
construction. In our work we adopted one of possible models of the S Scuti star 
XX Pyxidis calculated by Pamyatnykh et al. (1998). Our choice is motivated, 
in part, by the interesting data that we have for this object and, in part, by 
the fact that the star is a moderately evolved object, just like most of the well 
studied S Scuti stars. 



2Tn 



(69) 
(70) 
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The model mass is 1.9Mq, radius 2.O6i?0, effective temperature SOAbK, 
and luminosity IS.SLq. The model is slightly evolved, the central hydrogen 
abundance is 0.39. No rotation is assumed. 

We computed nonadiabatic acoustic modes in our model. The unstable 
modes of l<2 are found in the frequency range 3.7<c7<6.2. The most important 
parameters of the modes are listed in Table 1. Additionally, the growth rates 



Table 1 
Parameters of unstable acoustic modes 
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of the acoustic modes versus their frequencies are shown in Fig. 6. We will also 

use another linear characteristics of the modes, the complex factor /, which is 
the ratio of relative flux to relative radius variations. The absolute values of 
this factor are close to 30 for all the modes listed in Table 1. 

5.2 Gravity Modes 

The properties of g-modes are determined by the behavior of the Brunt- Vaisala 
frequency. It is presented in Fig. 7. Here we summarize the more general 
discussion presented by Dziembowski (1982) and Dziembowski and Krolikowska 
(1985). 

Gravity modes can propagate in a region where their frequencies are lower 
than the Brunt- Vaisala frequency (g-mode propagation zone) and this region 
have to be surrounded by evanescent zones where the Brunt- Vaisala frequency 
is lower. The frequencies of the g-modes that are efficiently coupled to the 
acoustic mode are close to one half of its frequency. Thus, as we may see in 
Table 1, the interesting gravity modes have frequencies between about 1.8 and 
3.1. We see in Fig. 7 that there are two separate propagation zones for such 
frequencies, one extending from r/R^O.l to r / and another in the 
range 0.95 <r/i?< 0.98*. Therefore, we have to consider two sets of g-modes of 

*In fact there is another, very narrow, propagation zone at r/R<0.1 for modes with iT^3, 
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Fig. 6. Dimensionless growth rates of the unstable acoustic modes of low-l in our model versus 
their dimensionless frequencies. 



vastly different properties. 

The inner g-modes can be studied in the asymptotic approximation (sec, 
e.g., Dziembowski 1982, Van Hoolst et al. 1998). It is valid when the wavelength 
of the mode is much shorter than the local hightscale. The radial wavenumber 
of a mode k is given by 



kk = — = -\U-k 
r r \ 



■'^BV 



(71) 



where Afsv is the dimensionless Brunt- Vaisala frequency, and Afe = Zfe(Zfc + l). 
The displacement, r?^., is expressed in the form 



where the eigenvector 

(r) = {yi,k {r)er + -2fc(r) V^) F™"' {0, ^) . 
The radial eigenfunctions in the asymptotic approximation are given by 

fk{r) 



yi,k{r) = 



fdfkir) 



.A-Va 



but this has no practical meaning. 



(72) 
(73) 

(74) 
(75) 
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Fig. 7. The dimcnsionlcss Brunt- Vaisala frequency versus radius in the adopted model. The 
atmospheric values are about 30. 



where 

fk{r) ='Dk\/ r/Kk sinipk, '^''^ j '^''^^^'^^ 
the structure functions A and Vg are given by 



V„ 



V _ 1 dlnp 
r ~~f dlnr' 



A 



dlnp 
dlnr 



(76) 



(77) 



r is the adiabatic exponent and remaining symbols have traditional meaning. 
When the mode degree Ik is high so is the radial wavenumber, implying high 
radial order, Uk- In such small change of the frequency is sufficient 

to change the order n by one, which means the spectral density of g-modes is 
very high. Thus, for each / it is quite easy to find a pair satisfying the resonant 
condition. We find in our model that the frequency distance between inner g- 
modes of a given I and consecutive n's, for the modes that are interesting for 
us, can be approximated by 



S(Tl,n = <^l,n-(ri,n+l'' 



0.36o-;,„ 

I ■ 



(78) 
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The inner g-modes are quasiadiabatic which means the only important nona- 
diabatic quantity is the damping rate. In this approach it may be approximated 

by 

-lk = -^, (79) 

where Ttu is the thermal timescale of the inner g-mode cavity. In our case it is 
roughly equal to 10^°, which implies the damping rates of the relevant modes 
range from 3 x 10~^ at 40 to 2 x 10~^ at 1000. Thus, the stability criterion 
given by Eq. (17) is not satisfied, because the driving rates of the acoustic 
modes, given in Table 1, are typically much higher. It is necessary to take 
into account at least several g-mode pairs to satisfy Eq. (48). As was shown 
in Subsection 2.2.1, this implies that the coupling to the inner g-modes cannot 
result in constant amplitude p-mode pulsation. 

Neither asymptotic nor quasiadiabatic approximation can be applied to 
modes trapped in the outer g-mode cavity. These outer g-modes still have 
relatively high degrees I but, because of the narrowness of the cavity, the radial 
orders n are low, typically 1 or 2. The nonadiabatic damping rates of these 
modes are typically between about 0.02 at (Ti_2^2.1 and about 0.5 at 0-1^2^;; 2.8. 
The strong damping of these modes might suggest they are not likely to be 
excited. However, the coiipling coefficients are also orders of magnitude larger 
than in the case of the inner g-modes, and the critical amplitudes given by 
Eq. (7) can be comparable to those for the latter. Neither in this case we expect 
a stationary pulsation because the second of the stability criteria (Eq. 18) is 
not fulfilled because he have <C 1 and 70 <C |7i,2|- Since the criteria given by 
Eq. (17) as well as Eq. (48) with N =1 are satisfied and the damping rates are 
much higher than the growth rates, a single outer g-mode pair is sufficient to 
halt the acoustic mode growth. The instability of the stationary solution in this 
case implies that the system reaches a limit cycle solution. 

In cooler 5 Scuti stars p-mode driving extends to lower frequencies, such 
that the resonant g-modes have their frequencies below the Mbv minimum. In 
this case the g-modes propagate through the whole radiative interior. Clearly 
the spectrum of such modes is very dense and damping is intermediate between 
that of the inner and the outer g-modes. 



5.3 The Coupling Coefficient 

Various authors give various expressions for the coupling coefficient. While the 
form of equations is basically always the same, the coupling coefficient may be 
expressed in various equivalent ways whose equivalence is not always apparent. 
Moreover, different formulae have different numerical properties and it is not an 
easy task to choose the best one. We restrict ourselves here to adiabatic formulae 
for two reasons. Firstly, there is no nonadiabatic formulae in literature available 
to us. Secondly, in the inner cavity the modes are indeed nearly adiabatic, 
yielding the adiabatic coupling coefficient to be an excellent approximation. 
This is not true for the outer g-modes. Therefore, we should treat the results 
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obtained for them as rough estimates. We will also discuss the effects of possible 
negligible coupling to these modes. 

The most detailed derivation of the adiabatic nonlinear oscillation equations 
and the formula for the coupling coefficient is given by Van Hoolst (1993, 1994a). 
This formula clearly shows symmetry properties of the coefficient, but it is not 
good for numerical calculations in the case of high-order modes because there 
arc many big terms canceling each other. Numerically much better formula of 
Dziembowski (1982) has a problem with treating the sharp peak of the Brunt- 
Vaisala frequency slightly below r = 0.1ii (see Fig. 7). The most appropriate 
for us seems to be the formula derived by Kumar and Goodman (1996). Trans- 
formed to our notation it takes the following form 



L^/rsdy, (80) 



^=87rG 

where H is the coupling coefficient introduced in Eqs. (1),(2). The integration 
is done over the stellar volume, and the third-order Lagrangian is given by 

+ ^^^^pdiyUdiyQ^ - ^P;ia(div^,)2 -p;,^^div^,div^o - 

(^^s^^div^o + ?gdivC,?^) - \{P;^;r.k + P^■,^■,r.kWg^g■ (81) 

In the last expression is the acoiistic mode eigenvector and the symbol 
denotes the sum of the gravity mode eigenvectors. The semicolons denote the 
covariant derivatives, upper indices denote contravariant vector components, 

Pa = -^a- Vp-rpdiv|„ 

is the Eulerian pressure perturbation of the mode a, is the gravitational 
potential, and 'I/q is its Eulerian perturbation of the acoustic mode. We make 
use of the Cowling approximation for the g-modes, ^''^ = ^'2 = 0. 
We introduce eigenfunction 

y2,k = ZkCcrl, 

where C = 3{r/R)^M/Mr and is the mass in a sphere of the radius r^. 
Additionally, for the acoustic mode we introduce eigenfunctions ys o,y4 o defined 
by 

— = gvifiYo- 



dr 



tNote that A/C is the square of the dimensionless Brunt- Vaisala frequency Nbv- 
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Then wc substitute cigenfunctions into Eq. (81), keep only terms with products 
of the functions of three different modes (only such terms contribute to the 
resonant coupling), and after very long and tedious algebra we obtain finally 



H 



where 



j ^^|2/i,iyi,2 yifl(^A{A + Ccjl-U-Vg)+2Ca^cj2{2-Vg)^ + 
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(82) 



Mr 



If the acoustic mode is radial the linear Eulerian perturbation of the gravita- 
tional potential can be obtained analytically and the function 2/2,0 has a different 
meaning than for the nonradial modes. It changes a bit the final formula for 
the coupling coefficient. It turns out, however, that if we put 2/3,0 = and 
2/4^0 = —C^ 2/1,0, Eq- (82) is valid in this case, too. 

In the asymptotic approximation we have to substitute Eqs. (74)-(76) into 
Eq. (82). The products of trigonometric functions of phases ijji,2 can be pre- 
sented as sums or differences of trigonometric functions of ipi+ip2 and ipi—ip2, 
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e.5-, sin?/'iCOS'(/j2 = l/2[sin('!/'i+?/'2) + sin('0i — ■02)]- Wc ignore terms varying with 
the rapidly changing phase Tpi + ip2 which nearly cancels out upon integration. 
We see that only modes with not too large difference of the radial orders n 
couple effectively because only then the phase difference changes slowly. The 
final formula for the coupling coefficient is obtained after keeping the dominant 
terms in Ai_2 
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(83) 



To obtain the expression for the coupling coefficient, C, we still need the asymp- 
totic formula for the moments of inertia. Using Eqs. (73)-(76) in the general 
expression, 



we get 



4 



1 



(84) 



1 

We adopt the normalization of spherical harmonics such that the angular inte- 
gral in Eq. (84) is equal 47r. Together with the standard normalization of radial 
eigenfunctions, yifi{R) = l, this gives us a simple interpretation of the amplitude 
^0 (see Eq. 72). Its absolute value, Bq, is then the surface averaged value of 
the 6R/R amplitude. Now, looking at Eqs. (3), (83), (84), we see that the coef- 
ficient C is independent of Pi, 2- It is also almost independent of gravity mode 
degrees, ^1,2, in the limit of ^1,2 ^ ^o- The dependence on Iq and mo remains. 
The latter comes only from the angular integral Z (see Dziembowski 1982 for 
explicit expressions). 

Only some sets of values of angular degrees give the nonzero coupling. The 
case of Z = takes place if Ik - Zi| > Zo or ii (-l)'^-'! ^ (-1)'°. Since we study 
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acoustic modes with lo<2 we need to consider the cases l2 = h for Iq = modes, 
^2=^1 + 1 for /o = l modes, and l2=li and l2=li+2 for Iq=2 modes. Additionally, 
the azimuthal numbers m have to satisfy the condition mo = mi+m2- 

In our numerical calculations for each acoustic mode we consider a set of 
inner g-mode pairs of h in the range from 40 to 1000, and for each li the range 
of An=n2—ni is from —20 to 20. Furthermore, for each (Zi,An) pair we consider 
various (mi,m2) pairs giving nonzero coupling coefficient. The dependence on 
li is indeed weak and is illustrated in Fig. 8 for the lowest unstable radial mode 
and lo = 2, mo = mode with cto = 3.883. 



0.003 I J, I I 0.015 I 




Fig. 8. The coiipliug coefficient C for the inner g-mode pairs and the lowest unstable radial 
mode (left panel) and io = 2, mo = 0, <to=3.883 mode (right panel). In the latter case /2 = h+2 
and the g-mode aaimuthal degrees are mi =m2 =0 which give maximum value of |Z| for these 
/-values. 



For the radial mode (left panel) most of the coupling comes from the region 
of a broad maximum of the Brunt- Vaisala frequency, 0.1 < r/i? 0.5 and if 
the modes 1,2 have different frequencies, their propagation zones have different 
sizes. These differences decrease with increasing li since the spectra densities 
increase. This is the cause of a bit different An-dependence of the coupling 
coefficients for different h. 

The two lowest lo = 2 modes (one of them is in the right panel) are in 
fact mixed modes with quite large values of eigenfunctions in the deepest part 
of the radiative region. Then, most of the coupling comes from the region 
r/i?< 0.1, and the coupling coefficient is almost independent of the upper limit 
of the propagation zone. Therefore, the coupling coefficient dependence on An 
is nearly the same for different h numbers. 

We can also see in the plot that if li = l2 (left panel) the coupling coefficient 
is symmetric about n2 — ni = 0. It is expected because the change of the sign 
of An means just interchange of the two g-modes and the coupling coefficient is 
symmetric with respect to it. Obviously, this is not true for li 7^/2 (right panel). 

The coupling coefficients have maxima near ni=n2 for each set of Iq,(Jq,12—Ii- 
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For the radial modes these maxima arc between 2.7 x lO'^ at ctq = 3.909 and 
1.3 X 10~^ at (To = 5.859. In the case of Iq = 1, mo = the maxima are between 
1.1x10-3 at 0-0=4.125 and 0.4 x 10-^ at a-o = 6.162. For Zo = l,|mi| = l the 
couphng coefRcicnts are about 30% larger due to larger values of Z. For Iq = 1, 
the values of coupling coefficients are a few times lower than in the case of radial 
modes due to lower p-mode amplitudes in the g-mode cavity. See also Fig. 9. 

As has already been mentioned, in the case of Iq = 2 we have to consider 
I2 = h together with l2 = li + 2. Typically, the l2 = h + 2 case is more important 
because the coupling coefficients are higher. The coupling coefficient maxima 
at mo = and mo ~ 1 arc between about 1.4 x 10^^ at a = 3.883 and 1 x lO""^ 
at (To = 5.793 while for mo = 2 the maxima are higher by a factor ~ 1.5. The 
coupling coefficients for the two lowest Zo = 2 modes are significantly larger than 
the coefficients for the remaining modes becaiise these two modes have mixed 
character and high values of eigenfunctions in the region around r/R = 0.1. 
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Fig. 9. The coupling cocfBcicnts maxima for the acoustic modes with m = 0. For Z = l,m = l 
the coupling cocfBcicnts arc about 30% higher than those for / = l,m = 0. For 1 = 2, m=l the 
coupling coefficients are similar to those at i = 2,m = 0. For Z = 2,m=2 the coupling coefBcients 
are higher about the factor ~1.5 than those for ( = 2,m = 0,l. 

In the case of the outer gravity modes we have to use Eq. (82) for H and the 
general expression for the moment of inertia. The angular degree dependence is 
the same as for the inner g-modes. It turns out that only modes with frequen- 
cies cri^2~o'o/2 and ni.2 = 1 are important. The estimated values of coupling 
coefficients for the modes in Table 1 range from about 7 at cto < 4 to about 
40 at CTo > 5.5. We should stress that these values should be treated as rough 
estimates, thus we do not give precise values. 
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5.4 Detuning Pciranieters and Critical Amplitudes 

Apart of damping rates and coupling coefficients we need values of the detuning 
parameters, Act. We notice first that in the absence of rotation the mode 
frequencies are independent of m. This means that the detuning parameters are 
independent of m, too. It is important because each set lQ,mQ,li,l2 characterizes 
whole multiplet of g-mode pairs. The pairs within the multiplet differ by m 
numbers of g- modes, but each pair has the same value of m\+m2- The pairs in 
the multiplet have the same values of Act, so, in view of what was said previously 
about synchronization, we may treat whole multiplet as a single pair. If 7^ 
the pairs within the multiplet have different coupling coefficients, so we take 
into account only pair with the maximum value of the coefficient. 

The inner g-modes with given angular degree, ^, and different radial orders, n, 
form dense frequency spectra. In such a case the determination of the detuning 
parameter would require an unrealisticly precise knowledge of the eigenfrequen- 
cies. Instead, we adopt statistical approach. For each g-modes' angular degree, 
I, we make a random choice of position of the frequency spectrum relative to 
the frequency of the acoustic mode. Then we compute detuning parameters 
for various g-mode pairs. The relatively strong coupling takes place only for 
ni «n2 so we may neglect |ni — n2| > 10 (see Fig.8). This means also that the 
relevant pairs have cti « CT2. Moreover, we explain below that only one pair for 
each {liil^) is important. 

At given I the frequencies of g-modes with consecutive n may be treated as 
equidistant in a relevant range of radial orders. Since li^l2 and nif«n2, we 
have from Eq. (78) Sui^^m ^6'^i2,n2i ^^nd 

ACTii,„i + l,(2,„2_l = CTo — (CTi^,„i+i -|-CT/2,„2-l) ~ 

~ CTq — (CTii,„i — <5CTj^,„i +Cr(2,„2 +(5CTj2,„2) PS 

K. Act„i,„2. (85) 

This implies that all the relevant pairs with even rii — 77-2 have close detuning 
parameters and may be treated as synchronized. As was explained in Section 4 
such pairs may be treated as one effective pair with the coupling coefficient 
equal to the maximum value among even ni — n2 pairs. Obviously, the same is 
true for odd n\ — n2 pairs. Finally, because 

|min(ACT)even m -ria " inin( ACT)odd ni -n2 I ~'^Crii, 2,^1,2 

is much higher than 70 for the relevant g-mode pairs, i.e., those with I of the 
order of a few hundreds (Eq. 78), we may take into account only this pair which 
has smaller |Act|. The other one has too high detuning parameter and is not 
excited. 

Even the lowest value of |Act| is often much higher than 70, and then we 
totally neglect g-mode pairs at a given {hjh)- What remains and is taken to 
further computations are single pairs at some values of h. 

It turns out that the most important pairs are in the range 100< Zi,2;$300. 
The lower-Z pairs are negligible for two reasons. First, their spectrum is sparse 
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and few of them have Act ^70. Second, they have very small damping rates, 
which implies they hardly contribute to the overall damping (Eq. 48). On the 
other hand the pairs with high Zi 2 are strongly damped and turn out to be 
inactive in time-dependent solutions. 

When the rotation is taken into account, the number of independent pairs 
increases due to splitting of g-mode frequencies. This effect will be studied in 
more details in Subsection 5.5.2. 

In the case of the outer g-modes the damping rates are much higher than the 
acoustic mode driving rates. This implies that only one pair suffices to halt the 
acoiistic mode growth and the system reaches the limit cycle solution, similar 
to that studied in Subsection 3.1. In this case the critical amplitude (Eq. 7) is a 
good estimate of the average acoustic mode amplitude. This critical amplitude 
is insensitive to the detuning parameter because it is much smaller than the 
damping rate and may be neglected in Eq. (7). This also implies that taking 
into account rotation does not change the expected acoustic mode amplitude. 

The outer g-mode pairs have to be taken into account if the average acoustic 
mode amplitude determined by the interaction with only the inner g-mode pairs 
is higher than the critical amplitude for the excitation of the outer g-mode pair. 
In particular, such a situation takes place for the acoustic modes with relatively 
high growth rates, because then many inner g-mode pairs are necessary to halt 
the driving. The detuning parameters of those pairs are also relatively high. 
Even though there is no simple dependence of resulting average acoiistic mode 
amplitude on detuning parameters of active inner g-mode pairs, the qualita- 
tive dependence is similar as in the single-pair case, i.e., the higher detuning 
parameters are, the higher the amplitude is. In our model, higher frequency 
acoustic modes have high growth rates and we expect high average amplitudes 
determined by the interaction with the inner g-modes. Then, the outer g-mode 
pairs may dominate the interaction. 

5.5 Time Evolution 

5.5.1 The Case of Slow Rotation 

In the case of sufficiently slow rotation, the splitting of the g-modes is smaller 
than their damping rates and the multiplets may be treated as single modes, as 
was explained in Subsection 3.3. Here we study a few examples of such systems. 
In all cases we set nio—O. Other values of mo give slightly different values of cou- 
pling coefficients, which results only in slight change of the amplitudes. Indeed, 
if we increase the coupling coefficients in Eqs. (30)-(32) by some factor and de- 
crease the amplitudes by the same factor, the equations will remain unchanged. 
This means the amplitude values scale with coupling coefficients without change 
of the time dependence which is governed by the p-mode growth rate. 

At first, we consider the = 2,cro = 3.883 mode together with an ensemble 
of several hundred inner g-mode pairs in the ?i-range between 40 and 1000. 
The beginning of the time evolution of the acoustic mode amplitude is shown 
in Fig. 10. The initial value of the amplitude is set to be equal to the low- 
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est critical amplitude. The initial g-mode energies were randomly distributed 
between zero and the acoustic mode initial energy. At the very beginning the 
acoustic mode amplitude drops as the result of the interaction with hundreds 
of artificially excited g-modes. Then, the modes become temporarily uncoupled 
and g-mode amplitudes fall down, while the acoustic mode amplitude starts 
to grow exponentially. Next, several g-mode pairs become excited due to the 
interaction with the acoustic mode and a statistical equilibrium is reached after 
about 5 X 10'* days. The average amplitude in this state is about 0.001 and 
average modulation period is about 3000 days. It is of the order of l/7o«5600 
days. Only inner g-modes actively interact with the acoustic mode. 




5x10* 105 

time [days] 

Fig. 10. The time evolution of the lo = 2,ao = 3.883 acoustic mode amplitude. 

The ^0 = 2,(70 <4 modes are not typical unstable modes in our model. They 
are, in fact, mixed modes with significant amount of energy concentrated in 
the internal g-mode cavity. This causes strong coupling to the inner g-modes 
and weak linear driving, the latter being the result of a large mode inertia. 
Therefore, we performed the same time integration for the radial ao = 3.909 
mode. Although the outer g-modes should be taken into account in this case, 
we first consider only inner g-modes in order to compare the results with the 
previous case when the growth rate was significantly smaller. 

The time evolution of the radial mode amplitude in the statistical equilib- 
rium without outer g-modes is presented in Fig. 11. We see that the average 
amplitude is about 0.01 and the typical modulation period is about 1000 days. 
This timescale is again of the order of l/7o~1600 days. The amplitude is about 
10 times higher than in the case in Fig. 10, mainly due to coupling coefficients 



33 



which are smaller by a factor of 5. The remaining factor of 2 turns out to be 
the result of dependence of the average amplitude on the growth rate. This 
dependence is expected because 70 determines Act range of active g-mode pairs 
and, as in the case of a single pair, there is some dependence of the average 
amplitude on detuning parameters. 




3.5x105 3.6x105 3.7x105 

time[days] 

Fig. 11. The time evolution of the lo = 0, ctq =3.909 acoustic mode amplitude without outer 
g-modes. 

The time evolution of the same radial mode with the outer g-modes taken 
into account is presented in Fig. 12. There is a qualitative difference in ampli- 
tude behavior between this case and the previous one. This difference is caused 
by the presence of outer g-modes which are strongly damped and strongly cou- 
pled to the radial mode. The sudden jumps of amplitude, lasting less than 10 
days, coincide with very short time intervals during which the outer g-mode is 
excited. 

Different behavior of the modes in this case than in the previous ones is 
caused by orders-of-magnitude larger ratio |7i,2|/7o of the dominant g-modes. 
Some of the inner g-mode pairs also remain excited, which is manifested by 
smooth amplitude changes, e.5., around 3.12 x 10^ days. However, the inner 
g-modes play a secondary role in this case. In spite of different behavior of 
the amplitude, the modulation timcscale is similar to that in the case without 
inner g-modes. This is the manifestation of the fact that also in this case the 
timescale is governed by the acoustic mode growth rate. 

This case resembles the case of a single g-mode pair, because only one pair 
is sufficient to stabilize the acoustic mode. The average amplitude of the latter 
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Fig. 12. The time evolution of the lo =0,ao = 3.909 acoustic mode amplitude with outer 
g-modes. 

is of the order of the critical amplitude to the excitation of the outer g-mode, 
equal approximately 0.005. This value is smaller than the average amplitude 
in the previous example, i.e., determined by inner g-modes and equal ~0.01. 
Therefore, the interaction with the outer g-mode dominates. 

The last example is the case of the Iq = 1,(To = 5.495 acoustic mode which 
is one of the most strongly unstable modes and for which the coupling to the 
inner g-modes is one of the weakest. The time evolution of this mode is shown 
in Fig. 13. Only an outer g-mode pair remain effectively interacting after several 
hundred days. The system converges to the periodic limit cycle with this pair. 
Very short modulation timescale, about 100 days, is caused by very high growth 
rate, l/7o~50 days. 

The high value of 70 demands many weakly damped inner g-mode pairs to 
stabilize the system. Moreover, these modes are very weakly coupled to the 
acoustic mode and, consequently, the resulting average amplitude of the latter 
would be much higher than ~0.02, i.e., the critical amplitude for the excitation 
of the outer g-mode pair. Therefore, the outer g-mode pair strongly dominates 
the interaction and we cannot even notice the presence of the inner g-mode 
pairs. 

Other unstable acoustic modes have average amplitudes similar to those 
in the last two cases, ie.,of the order of 0.01 -j- 0.02, and the interaction is 
dominated by the outer g-modes. 
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Fig. 13. The time evolution of the lo ■ 
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= 1, (TO = 5.495 acoustic mode amplitude. 



5.5.2 The Case of Moderate and Fast Rotation 

Since the modes of close frequencies effectively act as independent only if the Acr- 
separation is sufficiently large, rotational splitting may significantly change the 
picture only if rotation is fast enough. In this section we will try to determine 
the critical rotation rate at which modes can be treated as independent. In 
Subsection 3.3 we showed that the separation must exceed 7 for pairs to be 
regarded independent. 

Now, let us consider Act separation between members of g-modc multiplets 
that may be coupled to low-degree acoustic mode. Very good approximation 
for the high-Z inner g-mode's eigenfrequency in the corotating frame, up to 
quadratic terms in rotation is (Dziembowski et al. 1988) 



*^rot f^nrot 



i2_ ' 



(86) 



where flmi is dimensionless rotational frequency of a star. The g-modes in each 
pair have angular degrees h ~ I2 and mi ~ — TO2 so we neglect the difference in 
linear terms in rotation in the formula for the detuning parameter 



ACTr, 



' AcTnrot -I 

CTO 



4 / miflrot 



(87) 



where we used Cnrot ~cro/2. The rotational frequency change of the acoustic 
mode is unimportant here, since it only shifts all the detuning parameters by the 
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same value, without changing their statistical properties. The rotation begins 
to be important when 

^L>^- (88) 

In the numerical experiments from the previous subsection with only inner g- 
modes taken into account, the typical I of active pairs was 200. Thus, from 
Eq. (88) we find that the rotation becomes important at Orot ^0.001. Assuming 
uniform rotation of the star, this value in our model corresponds to equatorial 
velocity below 1 km/s. 

The linear terms in rotation are higher than the quadratic ones for very 
small rotation rate and h^h- However, at lo = 2 and /i~10^ we estimated that 
the linear Aa-splitting exceeds the quadratic one at rotation rate smaller than 
0.4 X 10~^. At such low values the splitting is smaller than |7/| so the rotation 
is negligible. When the rotation becomes significant, only quadratic terms play 
role. 

In order to assess terminal amplitude of acoustic mode in the presence of ro- 
tation we performed a number of numerical experiments. To reduce the number 
of g-mode pairs we divide each multiplet into groups of pairs such that detuning 
parameters in each group diff'er less than the value of the damping rate, 

|Ao-;,„-A£7i,™/|<|7;|. (89) 

Each group of pairs is treated as a single pair with mean detuning parameter 
and coupling coefficient equal the maximum one within the group. We ignored 
the influence of rotation on the coupling coefficients. The detuning parameters 
of different groups differed more than the damping rate of the multiplet. In 
particular, in the absence of rotation whole multiplet formed one group of pairs 
and was thus replaced by a single pair. 

An example of the time evolution of the Zq = 0,(To = 3.909 mode amplitude 
with the rotation frequency Vt^ot — -01 is presented in Fig. 14. We see that the 
average amplitude is about 0.0025. It is about factor of ~4 lower than in the 
case without rotation (see Fig. 11). Moreover, the characteristic timescale is 
somewhat longer than in the case without rotation. This is the result of much 
denser Act spectrum of active pairs. The dense spectrum implies that there are 
many independent low-Z pairs in the vicinity of exact resonance and there is no 
need to excite pairs with higher Act. Consequently, the modulation timescales 
are shorter, but still they are of the order of I/70. 

The transition from no rotation to fast rotation is shown in Fig. 15. At very 
slow rotation its influence on the average amplitude is negligible. The rotation 
becomes significant at f^rot ~ 0.002 which is of the order of the simple estimate 
obtained above. Further increase of the rotation rate causes that the average 
acoustic mode amplitude decreases. Above Jlrot ~ 0.01 the effect of rotation 
saturates. We may treat the rotation rate at which the average amplitude is 
reduced by the factor of 2 as the limit between slow and fast rotation. This limit 
in our case is Vl^ot ~ 0.003. Assuming uniform rotation this value corresponds 
to equatorial velocity about 2 km/s. 
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Fig. 14. The time evolution of the ao = 3.909 radial mode interacting with the g-mode pairs 
in the presence of rotation, Qrot =0.01. 
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Fig. 15. Average aiiiijlitude of the ao =3.909 radial mode interacting with the inner g-modes 
versus dimensionless stellar rotation rate. 
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If wc take into account outer g-modes it turns out that they dominate at 
slow rotation (see previous subsection). At fast rotation inner g-modes are able 
to limit amplitude growth of the radial mode at values below critical ones for 
outer g-modes which thus remain uncxcitcd. 

The role of rotation is important for the p-modes whose amplitudes are 
determined by the interactions with the inner g-modes. In our model they are 
weakly unstable ?o = 2,(To<4 modes. Their amplitudes are reduced by the factor 
of a few due to fast rotation. 

On the other hand, for high frequency p-modes (ctq > 4) even fast rotation 
cannot change the situation, because even at the fastest rate the outer g-modes 
must be excited and they determine the p-mode amplitudes. As was explained 
in Subsection 5.4, the outer g-modes are insensitive to rotation. Thus the values 
of the amplitudes of the (Tq > 4 p-modes are the same as in Subsection 5.5.1, 
i.e., of the order of 0.01 0.02. 

5.6 Luminosity Amplitudes 

In order to obtain the expected luminosity amplitudes, wo have to multiply ra- 
dius variations amplitudes by a nonadiabatic coefficient /. As was mentioned in 
Subsection 5.1, the absolute value of this coefficient for all considered acoustic 
modes is about 30. The effects of flux averaging over the stellar disc play role 
mainly for 1 = 2 modes, and they reduce effective |/| to the value of about 20. 
Moreover, for the modes whose radius amplitudes are determined by the inter- 
action with the inner g-modes, we adopt the values obtained at fast rotation. 

In the case of radial acoustic modes the evaluated bolometric luminosity 
amplitudes range from about 0.07 mag at cr = 3.909 to more than 0.5 mag at 
(7 = 5.859. Only the (7 = 3.909 radial mode has amplitude determined by the 
interaction with the inner g-modes and the value of the amplitude is calculated 
with the rotation taken into account. 

All the I = 1 acoustic modes have amplitudes determined by the outer g- 
modes. The values of the amplitudes range from ^0.15 mag at (t = 4.125 to 
more than 0.5 mag at a > 5. 

The 1=2, acoustic modes with (j<4 have relatively small expected luminosity 
amplitudes, about 5 mmag. This is the combined result of their low growth 
rates, strong coupling to the inner g-modes, fast rotation and disc-averaging 
effect. However, at tr > 4 the expected amplitudes are again high, of the order 
of a few tenths of a magnitude. 

Comparing obtained luminosity amplitudes with those observed in XX Pyx 
or other 5-Scuti stars, we see that only low-frequency I = 2 modes in our model 
have expected amplitudes below 10 mmag. Higher frequency acoustic modes 
have amplitudes equal a large fraction of a magnitude. These values are typical 
for HADS stars. Neglecting outer g-modes results in even higher amplitudes. 
This suggests that nonresonant saturation effects play an important role in limit- 
ing the amplitudes. However, at present stage we cannot exclude the possibility 
that our treatment of the outer g-modes is too crude and more realistic modeling 
of those modes would give better consistency with observations. 
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5.7 Effects of Saturation of The Linear Driving 

The saturation of the instabihty mechanism is manifested by the nonHnear de- 
pendence of the growth rates on mode amphtudes. The nonhnear growth rates 
may be written in the form (see, e.5.,Buchler and Goupil 1984) 



where Skj < arc saturation coefRcients and the summation is made over all 
modes excited in the star. 

The simplest way to consider the saturation mechanism is to reduce the 
values of the growth rates. Then, as implied by Eq. (48), the number of g- 
modes necessary to achieve the statistical equilibrium is reduced, too. The 
lower driving rate, as we explained in Section 4, means also that the Act range 
of active pairs is smaller. Consequently, the expected amplitude of the acoustic 
mode becomes lower and the modulation timescale is longer. 

The effect of the growth rate saturation is relevant only in the interaction 
with the inner g- modes. For the outer g- modes there is no reduction of the 
amplitude because the damping rate is always much higher than the growth 
rate. 

The cases of very small growth rates were studied by Dziembowski and 
Krolikowska (1985) and Dziembowski et al. (1988). The amplitude of the acous- 
tic mode is then given by the lowest critical amplitude (Eq. 7) for the whole 
ensemble of the resonant g-mode pairs. It is so, because the pair corresponding 
to this minimum is able to halt the acoustic mode growth. The system reaches 
the stable equilibrium, described in Section 2.1.2, or a periodic limit cycle such 
as described in Section 3.1. The expected amplitude of the acoustic mode in 
both situations is close to the critical one. 

For the XX Pyx the simple assessment based on Eq. 7 with the rotation 
and strong saturation taken into account and the respective / coefficient leads 
to amplitude ranging from about 1 mmag at I = 2,a < 4 modes up to about 
0.1 mag at Z = 1,ct > 5 modes. The amplitudes of the high frequency acoustic 
modes arc still much larger than observed in XX Pyx. The fact that we find 
amplitudes considerably higher than Dziembowski and Krolikowska (1985) did, 
is a consequence of evolved 5 Scuti star model while they used a ZAMS star 
model. In our case the inner g-modes are more concentrated at the peak of the 
Brunt- Vaisala frequency below r = 0.1i? and consequently they are much more 
weakly coupled to the p-mode. The interaction with the outer g-modes is also 
weaker in our case because of higher damping rates. 

The unavoidable conclusion is that the dominant effect responsible for the 
amplitude limitation in 5 Scuti stars must be nonlinear saturation. In such 
a situation our basic approximation that acoustic modes remain uncoupled is 
not allowed. The coupling between various acoustic modes occurs through the 
collective saturation of the driving mechanism, as described by Eq. (90). The 
decrease of amplitudes of certain modes gives the chance for other modes to be 




(90) 
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excited. The net effect depends on the linear growth rates and the saturation 
coefRcients. In Fig. 6. we see that the hnear growth rates initially grow rapidly 
with frequency. However, the saturation coefScients are also growing with fre- 
quencies as both observations and hydrodynamical simulations show. The high 
frequency modes are both easy to excite and easy to saturate. It requires nu- 
merical calculations to predict the outcome of the collective saturation. The 
only thing that can be said at this stage is that the amplitude changes should 
be faster for high frequency modes. 

6 Conclusions 

We studied the resonant coupling between unstable acoustic modes and pairs 
of stable g- modes. We reminded the general analysis of the interaction given 
by Dziembowski (1982) and Dziembowski and Krolikowska (1985). The most 
important properties of the interaction are such that a g-mode pair gets excited 
if the acoustic mode amplitude exceeds certain critical value, and that a three- 
mode equilibrium solution is stable if the acoustic mode driving is weaker than 
g-mode damping and if the detuning parameter is not too small. We also showed 
that the nonadiabacity of the nonlinear coupling does not change the problem 
qualitatively. 

In the case of many g-mode pairs the equilibrium solution involving more 
than two pairs in general does not exist. Instead, we have a kind of statistical 
equilibrium involving many pairs exchanging energy with the acoustic mode. 
We showed that such a state may be stable if the sum of damping rates of 
active g-mode pairs is higher than the acoustic mode growth rate. Moreover, 
these g-mode pairs should be separated in the space of dctiming parameters by 
more than their damping rates. The pairs that have close detuning parameters 
get synchronized and interact effectively as a single pair. 

The separation of detuning parameters and the condition of the stability 
imply that the range of the detuning parameters of active pairs is of the order 
of the p-mode growth rate. The inverse of the growth rate determines the 
timescalc of the modulation of the acoustic mode amplitude. The fast amplitude 
variations cause that the expression for the critical amplitude derived in quasi- 
static approximation is no longer applicable and many pairs are excited despite 
being stable according to the parametric instability criterion applied to the 
time-averaged p-mode amplitude. 

We applied our theory to a seismic model of a ^ Scuti star XX Pyx. We 
consider all unstable acoustic modes of Z < 2. Their dimensionless frequencies, 
v^'k/[G<p>), are between, 3.7 and 6.2, corresponding to radial orders n>3. 
Such modes couple to g-mode pairs that propagate either in the outer layers, 
0.95 < r/i? < 0.98, or in the interior, 0.1 < r/i? < 0.5. The two groups of g- 
modes have very different properties, in particular the frequency spectra. While 
the number of the inner g-modes that can be excited due to the parametric 
instability is large, typically only one outer g-mode pair has to be taken into 
account. The outer g-modes are orders of magnitude more strongly damped 
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but also much more strongly coupled to the acoustic modes, resulting in critical 
amplitudes comparable to those for the inner g-modes. 

The 1 = 2 acoustic modes of low frequencies, o-<4, reach statistical equilibria 
determined by the interactions with the inner g-modes of 10^. The average 
amplitudes are below the critical values for the excitation of the outer g-modes 
which thus remain unexcited. The rotation causes that the average amplitudes 
are by the factor of a few lower than in the case without rotation. It is so because 
the rotation causes splitting of detuning parameters and above certain rate 
the pairs within multiplets loose synchronization and act independently. This 
increases the number of effectively interacting g-modc pairs and leads to lower 
acoustic mode amplitudes. The effect of rotation becomes significant already 
at equatorial velocity of the order of 1 km/s. Then, the average bolometric 
luminosity amplitudes are of the order of a few milimagnitudes. The typical 
modulation timescales are of the order of 10"^ 10'* days. 

The 1 = 0, cr = 3.909 mode has average amplitude determined mainly by the 
interaction with the outer g-mode pair, but the time behavior is affected by the 
presence of a few excited inner g-mode pairs. Effect of rotation is to increase 
the number of inner g-modes and these modes take control of the acoustic mode 
evolution. The average amplitude is now lower leading the outer g-modes unex- 
cited. The rotation again turns out to be important already at rates ~ 1 km/s. 
Then the average bolometric luminosity amplitude is about 0.07 mag and the 
modulation timescale is of the same order as in the previous case, i.e., 10"^ days 

The amplitudes of high frequency acoustic modes, cr > 4, are determined 
solely by the interaction with the outer g-mode pairs. Typically, the system 
reaches periodic limit cycle solutions. The rotation has no effect in this case. 
The predicted bolometric luminosity amplitudes are large, much higher than 
observed. The modulation timescales are of the order of 10^ -t- 10^ days. 

We found that the resonant mode coupling alone cannot be the main effect 
responsible for amplitude limitation in 6 Scuti stars, at least in hotter and 
evolved ones as XX Pyx. Our calculated luminosity amplitudes are too high, 
in particular for high frequency modes. There is an obvious conflict with the 
observations. 

A nonlinear saturation of the linear driving effect leads to a reduction of the 
acoiistic mode amplitudes. If strong enough it may lead to constant amplitude 
pulsation such as discussed by Dziembowski and Krolikowska (1985). However, 
in our application to an evolved 6 Scuti star model the calculated amplitudes 
are still higher than observed. The difference between the two results is a 
consequence of the difference in the internal structure of the adopted models. 
We stressed that in the case when the dominant amplitude limitation mechanism 
is the nonlinear saturation, the coupling between unstable acoustic modes must 
be taken into account. 
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